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Abstract 

A new quantum approach is presented that can account for the 
description of small ferromagnetic particles magnetization tunneling. 
An estimate of the saturation value of an external applied magnetic 
field along the easy axis is obtained. An analytic expression for the 
tunneling factor in the absence of an external magnetic field is deduced 
from the present approach that also allows to obtain the crossover 
temperature characterizing the regime where tunneling is dominated 
by the quantum effects. 
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Single-domain magnetic particles constitute a scenario in which macro- 
scopic quantum tunneling can be studied. In such systems of mesoscopic size 
the eletronic spins can form an aligned magnetic state that can be oriented in 
several directions. Quantum mechanics can then be called upon to estimate 
the tunneling probability that may occur, associated with the possibility of a 
change of the magnetization from one direction to another. The existence of 
an energy barrier to be tunneled that do not completely disappear even with 
a decrease in the temperature towards absolute zero has therefore been sug- 
gested and several attempts have been presented in order to estimate the rate 
of tunneling in these cases [TJ (2J [3l H] . In some approaches an adapted version 
of the WKB is proposed, while, from another perspective, other approaches 
make use of Feynman's path integral with su(2) coherent states from which 
the instanton method follows to treat this same quantum problem for a single 
spin in predicting the energy levels splitting . 

In general, those approaches start with a Hamiltonian that is associated 
with a model of the physical system and some considerations are presented in 
order to get a pair of complementary operators which should characterize the 
degree of freedom. To this end, the spin system which is characterized by a 
finite-dimensional state space, is converted into a continuous one by means of 
a transformation which leads to a pair of canonical operators satisfying a new 
commutation relation which is now dependent on the particle spin quantum 
number [5]. It is to be observed that this construction is valid for large spin 
values and that corrections of quantum character must be introduced for 
small/medium values. 

In the present paper we intend to show that a new quantum description 
recently proposed [6] can be used in such a way that the energy levels splitting 
in ferromagnetic particles can be obtained in a direct and simple way. 

We start from a model Hamiltonian describing the ferromagnetic particle 

® 

H = -KrSl + K 2 S 2 y + K 3 H Z S Z (1) 

written in terms of spin operators S y and S z with Ki, K 2 > 0, corresponding 
to the case of easy-axis anisotropy along the z axis and an external magnetic 
field pointing along the z direction. Putting this Hamiltonian in a form that 
emphasizes the creation and anihilation operators we get 

H = AS 2 +B(S 2 + + S 2 ) + CS 2 + K 3 H Z S Z , (2) 

where A = - (K l + K 2 /2) , B = -K 2 /A, and C = K 2 /2. From a formal 
point of view, it is interesting to see that this Hamiltonian is akin to that 
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one suggested to describe the Fe$ magnetic cluster (also with an external 
magnetic field) and that has been treated elsewhere [B]. 

The starting point of the proposed approach is the introduction of a 
quantum phenomenological Hamiltonian describing the spin system, written 
in terms of angular momentum operators obeying the standard commutation 
relations, and that reflects the internal symmetries of the system. It may also 
contain terms taking into account external applied magnetic fields. The de- 
gree of freedom that undergoes tunneling is considered a particular collective 
manifestation of the system, and it is assumed to be the only relevant one. At 
the same time, the temperature of the system is assumed so conveniently low 
that possible related termally assisted processes are not taken into account 
so that only quantum effects are considered. 

The next step consists in getting a new Hamiltonian that is an approxi- 
mate version of Eq. ([2]). This new Hamiltonian is obtained through a series 
of transformations performed on the matrix generated by calculating the ex- 
pectation values of Eq. (j2J) with the su(2) coherent states \j, z), where z is a 
complex variable and j characterizes the angular momentum state multiplet 
[7]. The Hamiltonian in the overcomplete spin coherent states representation 
is then given by 



the also known generator coordinate energy kernel j8]; it embodies all the 
quantum information related to the system we want to study. The procedure 
of extracting a new Hamiltonian - written now in terms of an angle variable 
- from Eq. ([3]) has been already shown elsewhere and will not be repeated 
here [6j [9l [10]. In fact, as it was proved there, the variational generator 
coordinate method can in this case be used to rewrite the Hamiltonian from 
which we start in an exact and discrete representation which can then be 
conveniently treated in order to give an approximate Hamiltonian in the angle 
representation. It is important to point that the Hamiltonian we obtain is 
in fact an approximate one, but we have also shown, by studying the Lipkin 
model [TTJ, that it is already a reliable Hamiltonian for spin systems with 
j > 5 [BJ. Furthermore, since this approach is based on quantum grounds 
from the beginning, it does not need to go through any quantization process. 
Also, it is not necessary to convert the discrete spin system into a continuous 
one as it is usually done [5j , at the same time that the quantum character of 
the angle-angular momentum pair is properly taken into account. 

Now, carrying out all the steps of the approach, we obtain the approxi- 




(3) 
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mate Hamiltonian in the angle representation, namely 



H{4>) 



1 d 1 d 

2d0M(0)d0 



(4) 



where 



M{<j>) 



1 



(5) 



2 (K 2 + cos 2 0) - ^H z cos 



is the effective mass, 



V (0) = -K l3 (j + 1) cos 2 + K 3 H z y/j(j + l)cos(f> 



(6) 



is the potential respectively, and j is the total spin of the ferromagnetic 
particle. By a direct inspection we see that the potential function has minima 
at = 0, 7r, while the maxima occur at 



respectively. It is immediate to see that for H z = 0, the maxima are at 
= vr/2 and 3tt/2. 

As has been mentioned before, the accuracy of this approximation can be 
checked by comparing the eigenvalues of the Schrodinger equation 



with those obtained by diagonalizing Hamiltonian (TjQ) in the {\jm >} state 
space. By its turn, it is immediate to see that the numerical results of Eq. 
(j7j) can be directly obtained from a Fourier analysis that give the energy 
eigenstates and the corresponding eigenfunctions. It has also been shown 
elsewhere [6] that already for j ' > 5 solutions of Hamiltonians of this kind 
are in good agreement with the ones coming from ([!]), being the deviation of 
the order of 1%. 

An interesting result can be immediately drawn from a direct analysis of 
the potential energy expression, Eq. fl§]). It is simple to see that the extrema 
of the potential energy depend on the strength of the external magnetic 
field. Furthermore, if we realise that an increase of that field tends to gather 
together the two maxima, at the same time that it rises the value of the 
minimum at = ir (or at = if the field is reversed) , we can expect to find 
the intensity of the field at which those three extrema colapse into a single 



K 3 H Z 




#(0) V k (0)=£A(0) 



(7) 
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maximum. At this particular value of the field there is a single minimum 
at = (or at = 7r if the field is reversed) and tunneling cannot occur. 
Therefore, using the proposed approach we can estimate the value of the 
external magnetic field at which tunneling no longer has any meaning if we 
consider the expression for the potential energy expression, Eq. ([6]), and look 
for the value of H z for which the minimum turns into a maximum. After a 
direct calculation we obtain 

so that 

#' im Kt r T 



is the value of the external paralel magnetic field whose multiples produce a 
matching of the energy levels in the two potential wells. The result obtained 
from a numerical diagonalization of Eq. ([1]) for the matching value in a test 
case is in good agreement with the predicted value of Eq. (jSJ). It is then 
evident that for an applied external field of magnitude if' im the resulting 
potential function presents only a single minimum. 

It is clear that if we want to discuss all the features of the energy levels 
associated with Hamiltonian (j3J) we may numerically solve the Schrodinger 
equation (E]) and get the results. Although this goal is simple and direct, it 
will not be pursued here. 

In what concerns the energy levels splitting associated with the tunneling 
through the potential barrier present in V (0) in the absence of the external 
magnetic field, i.e., H z = 0, we may take advantage of the even symmetry 
of V (0) as well as of M (0) to use the well-known WKB expression (h = 1) 



AS =— exp 

7T 



-Vzm 1 3 Vv (0) - £ 



(9) 



where u;& is the frequency associated with one of the minima of V (0), and 
M. is a constant suitably taken as an average of M (0) in the barrier region. 
Considering 

MU \ ^ d 2 V (0) 
M (0 min ) u b = — 



and taking into account that min = 0, ir and max = 7r/2, 3tt/2, we see that 

^ = 2y/Kl (Kl + K 2 )j(j + 1). (10) 



5 



On the other hand, the exponential factor can be treated in an analytic way. 
Taking 



V / -^ij(j + l)cos 2 0+ \S\ = VW\\1 



Kij (j + 1) 
\£\ 



■ cos z 



and calling k 2 = Kij (j + 1) / \£\ we see that the turning points of the integral 
are given by the equation cos (f) ex t = ±A; _1 . Therefore, 



1(E) 



cos^ 



fc2 - 1 F ( ■ fir \ 1 

— t I arcsm [k cos) , — 



k 



kE ( arcsin (k cos) , — 



if k > 1, 



where F (ip,p) and E (<f,p) are the elliptic functions of the first and second 
kind respectively [IB] . Since these functions are odd under the change <p — > 
— <p, and using the expression for cos0j jS , we obtain 



/(£) 



kEi k 



;ii) 



where K (|) = F (|, |j and E (|j = (|, |) are the complete elliptic 
integrals [IB] . Now, using their corresponding series representations we end 
up with the result 



7T /l 1 3 



2 8^ 3 64P 

so that, the integral to the lowest order reads 



(12) 



1(E) 



7T 7T 



\s\ 



2k 2\K l3 ( 3 + lY 
With that result we can then write the barrier penetrability factor 



(13) 



cxp 



2 1 1 V i^j (j + 1) 



(14) 



6 



In order to have an estimate of the tunneling factor we have to introduce 
the approximate constant mass. If we consider for instance the crude simple 
arithmetic mean in order to get an analytic expression 



M 



M(0 = O) + M(0 = f) 
2 



AK 2 (K x + K 2 ) 



Ki + 2K 2 



(15) 



for the constant mass, we obtain 



P ~ 



exp 




2K 1 j(j + l)K 2 (K 1 + K 2 ) 



[K x + 2K 2 ) 



On the other hand, the frequency expression can be directly obtained from 
Eq. ffTOl) so that the expression for the energy splitting, up to first order in 
the exponential, is 



As it can be immediately seen, P {£) = 1 at the top of the energy barrier 
(i.e., £ = Vmax = 0.0), thus we verify that this expression is not adequate to 
describe tunneling occurring in energy levels near or at the barrier top. In 
these cases, we may consider the Kemble-Hill- Wheeler/Miller-Good expres- 
sion [TU [151 [HH [IT] for the barrier penetration factor, namely 



where u t is the frequency associated with the concavity at the potential 
barrier top, and V = V ((ft = 4>top)- It is important to see that at £ = V we 
get P (Vq) = 0.5 instead of the usual WKB value P = 1.0, thus assuring that 
the Kemble-Hill- Wheeler/Miller-Good expression is a better approximation 
for this particular kind of situation. 

If we realise that the exponential factor is the dominant part of the energy 
splitting expression, Eq. we are lead to consider again expression <HM 
in order to estimate the temperature below which the quantum effects dom- 
inate over the thermal ones. Since the thermal ativation is proportional to 
exp(— Vq//cbT), where Vq is the barrier height, and in our case it is given by 




2K 1 j(j + l)K 2 (K 1 + K 2 ) 

(16) 



{K x + 2K 2 ) 



P{£) 



1 



1 + exp [f (V - £)] ' 
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Vo = Kij (j + 1), and considering that \£\ — Kij (j + 1) for j ^> 1, then the 
crossover temperature between the two regimes can be estimated by equating 

7T / 2M _ 1 

that gives 

J_ / K l3 (j + 1) 

The effective mass plays an important role in the above expression, telling us 
how the transverse anisotropy governs the behaviour of the relevant physical 
quantities related to the tunneling. In order to have an estimate of this 
behaviour let us consider Eq. f|T5|) for the effective mass. It is direct to see 
that for K 2 — * 0, M. —* oo so that T c — > 0, and consequently AS — > 0, as 
expected, since then there is no tunneling at all. It is to be stressed that 
the above quantitative results do not follow due to a particular form of the 
effective mass used. They would also follow for any reasonable expression 
for the mass. It is important to observe that the effective mass is the only 
term carrying the information related to the transverse anisotropy which is 
responsible for the quantum tunneling effects. 
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